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LONG-TERM  GOAL 

To  use  an  iterative  time-reveral  technique  to  provide  robust  focusing  along  the  ocean 
bottom,  with  little  signal  processing  effort  involved  and  a-priori  information  on  the 
environment,  to  enhance  the  detection  and  localization  of  a  proud  or  buried  target  in 
shallow  water. 

OBJECTIVES 

The  objective  of  the  proposed  MPL  component  of  a  SWAMSI  continuation  will  develop 
bi-static  applications  of  the  iterative  time  reversal  techniques  to  focus  along  the  ocean 
bottom.  Iterative  time-reversal  provides  a  simple  solution  for  self-adaptive  focusing  on 
strong  reflectors  (i.e.  scattering  targets)  located  on  the  ocean  bottom  without  relying  on 
predictive  or  modeling  capabilities  of  the  environment  and  of  the  target  of  interest.  These 
robust  focusing  properties  are  crucial  for  mapping  large  and  uncharted  areas  with  little 
signal  processing  effort  involved.  With  the  experiments  planned  by  Schmidt  et  al,  there 
will  be  enough  data  to  understand  the  ultimate  limitations  of  the  proposed  self  adaptive 
methods. 

BACKGROUND  FOR  FOCUSING  CONCEPT 

The  robust  focusing  and  pulse  compression  provided  by  time  reversal  techniques  can  be 
exploited  to  enhance  active  sonar  processing  in  shallow  water  environments.  A  successful 
time-reversing  array  (TRA)  retrofocus  can  be  created  without  any  prior  knowledge  of  the 
environment,  the  probe  source  or  the  reflector  location,  or  the  array's  transducer 
locations.  A  TRA  may  be  of  nearly  any  size  or  shape  and  may  operate  in  any  frequency 
range.  Moreover,  TRAs  with  sufficient  aperture  perform  well  in  multipath  environments. 
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Figure  1:  Diagram  of  an  Iterative  Time-Reversal  process  in  a  classical 
monostatic  configuration.  1  iteration  is  used  to  enhance  focusing  on  a  reflector 
located  above  sediments,  (a)  Broadside  illumination  with  a  pulse  signal,  (b) 
Backscattered  signals  are  recorded,  (c)  A  time  window  containing  the  potential 
target  signal  is  selected  in  order  to  be  time-reversed,  (d)  The  resulting  wave 
illuminates  the  bottom,  (e)  The  backscattered  signals  reveal  that  the  waves 
focused  on  the  target  since  the  target  echo  is  enhanced.  Further  iterations  can 
help  detect  weaker  targets. 


Thus  the  (iterative)  time  reversal  process  can  provide  a  means  for  self-adaptive 
focusing  in  the  vicinity  of  the  water-bottom  interface  in  order  to  improve  the  detection  of 
targets  laying  on  the  seafloor  or  buried  in  the  sediments  (see  Figure  1).  By  applying  this 
time-reversal  process  iteratively,  we  construct  a  set  of  signals  for  the  TRA  which  focus 
on  the  strongest  scatterer  located  in  the  selected  range  interval.  This  is  achieved  by 
repeating  the  iterative  time-reversal  process  over  a  sliding  time-gated  window  portion  of 
the  initial  reverberation  signals  recorded  in  a  waveguide,  thus  allowing  us  to  scan  the 
bottom  surface.  Reflectivity  maps,  similar  to  a  sonar  map,  but  with  an  enhanced  contrast 
for  the  strongest  reflectors  (or  scatterers),  at  the  water-bottom  interface  can  be 
constructed  using  a  simple  TRA.  When  the  reflectivity  contrast  between  the  buried  target 
and  the  bottom  reverberation  is  sufficient,  focusing  energy  on  the  target  is  accomplished 
with  very  little  signal  processing  effort.  Otherwise,  an  eigenvector/eigenvalue  analysis  of 
the  backscattered  signals  can  be  applied  to  isolate  the  target  contribution  and  to  construct 
appropriate  time-reversed  signals  that  focus  on  the  target  solely.  Hence  the  iterative  time- 
reversal  process  should  provide  a  simple  solution  for  separating  the  echoes  of  the  targets 
from  the  reverberation  of  the  ocean  bottom. 
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In  support  of  our  potential  iterative  and  ultimately  bi- static  concept,  we  mention 
the  previous  passive  iterative  time  reversal  implementation.  Only  the  knowledge  of  the 
inter-element  impulse  response  matrix  for  TRA  is  required,  i.e.  measuring  the 
backscattered  Green’s  functions  (or  impulse  response)  between  all  pairs  of  elements  of 
the  TRA  after  round-trip  propagation  between  the  TRA  elements  and  the  ocean  bottom, 
thus  including  the  whole  bottom  reverberation  response.  The  acquisition  of  the  inter¬ 
element  impulse  response  matrix  is  the  only  active  part  of  the  processing.  The  full 
physical  propagation  characteristics  between  the  TRA  and  the  ocean  floor  are  then 
contained  in  the  inter-element  impulse  responses  at  once,  and  no  extra  information  is 
required.  The  various  steps  of  the  time-reversal  process,  as  described  in  Fig.  1  can  then  be 
implemented  passively  by  simply  performing  correlations  of  these  inter-element  signals 
with  the  time-gated  window  of  the  reverberation  signals  selected  at  the  previous 
iterations.  This  allows  us  to  construct  passively  the  reverberation  time-series  of  each 
TRA  element  as  if  successive  back  and  forth  propagation  between  the  TRA  and  the  ocean 
floor  had  occurred.  Then,  the  time  reversal  iterations  can  be  implemented  passively. 

This  passive  iterative  time-reversal  technique  was  applied  to  detect  a  target 
(composed  of  8  small  glass  spheres)  with  a  curtain  source  receive  array  (SRA)  used  as  a 
time-reversal  array  (having  a  small  vertical  aperture  of  8  by  1.5m,  operating  from  3  to  4 
kHz)  that  is  a  surrogate  for  existing  sonar  systems  in  a  mono-static  configuration  (see  Fig 

1).  The  curtain  SRA  was  deployed  in  50m  shallow  water  and  the  target  was  located 
approximately  on  the  bottom  203m  away  from  the  SRA,  and  the  inter-element  response 
was  measured  over  3min.  This  reverberation  focusing  technique  appears  robust  with 
respect  to  environmental  fluctuations.  A  reflectivity  map  was  constructed  for  one 
iteration  of  the  time-reversal  process  (see  Figure  2. a  and  2.b).  The  target  appears  at  the 
estimated  range  as  a  bright  spot  in  the  main  diagonal  delimited  by  the  extent  of  time¬ 
gated  window  used  to  construct  time-reversed  reverberation  signals  (see  Figure  1). 
Detection  of  the  target  is  enhanced  over  a  conventional  broadside  return  (see  Figure  2.c). 

PROPOSED  RESEARCH  ITEMS 

1)  Develop  and  validate  synthetic,  iterative  time  reversal  aperture  for  SWAMSI 

2)  Relate  and  understand  the  performance  of  the  above  procedure  to  the  near  field 
measurements  that  are  being  performed  simultaneously 

3)  Extend  the  processing  method  to  a  cooperating  multi-platform  scenario 
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Figure  2:  a)  and  b)  Reflectivity  maps  for  a  target,  located  203m  away  from  the  SRA. 
These  maps  were  computed  over  sliding  time-gated  banning  window  portion  of  A=20ms 
of  the  reverberation  signals  (delimited  by  the  dashed  lines  along  the  main  diagonal)  for: 
a)  the  time-reversal  process  or  b)  1  iteration,  of  the  time-reversal  process.  The  horizontal 
axis  is  the  recording  time  (i.e.  the  round  trip  propagation)  converted  to  the  SRA/bottom 
distance.  d=t  co/2  (co=1515m/s)  between  the  SRA  to  the  bottom.  The  vertical  axis  is  the 
average  intentended  focusing  distance  along  the  bottom  being  scanned  by  the  TRA 
correponding  to  the  selected  time-gated  window.  Note  the  contrast  improvement 
achieved  for  the  strongest  reflectors  by  using  one  iteration.  The  target  location  is  the 
brightest  spot  at  203m.  c)  Comparison  of  the  intial  boradside  return  to  the  diagonal  of  the 
reflective  map  after  iterative  time-reversal. 


APPROACH 

NURC  MCM’06  Experiment:  Application  to  bi-static  detection  enhancement. 

Our  method  requires  the  addition  of  a  moving  (suspended,  drifting)  receiver  that  has  a 
common  sync  to  the  TOPAS  source.  The  previous  Passive  Iterative  Time  Reversal 
technique  can  be  applied  to  a  dynamic  bi-static  geometry  where  the  source  and  receiver 
roles  are  separate.  The  active  source  is  provided  by  a  TOPAS  parametric  source,  mounted 
on  a  telescopic  tower,  insonifying  a  target,  as  per  Schmidt  et  al.  At  the  same  time,  a 
moving  receiver  (e.g.  towed  hydrophone  from  a  surface  platform)  collects  the  bi- static 
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backscatter  signals  from  the  target  at  several  positions  along  its  trajectory.  These  data  are 
then  used  to  construct  a  synthetic  aperture  Time-Reversal  Array  which  can  be  used  for 
the  passive  iterative  time-reversal  process. 


Receiver  Trajectory 


Topas 

Tower 


Target 


Figure  3:  Schematic  of  bi-static  implementation  of  time-reversal  process  for  detection 
and  localization  of  Target  in  shallow  water.  The  Topas  tower  provides  successive 
illumination  of  the  target.  A  moving  receiver  (e.g.  towed  hydrophone  from  a  surface 
platform,  AUV...)  collects  the  bi-static  backscatter  signals  from  the  target  at  several 
positions  along  its  trajectory.  This  moving  receiver  can  then  be  used  to  construct  a 
synthetic  aperture  Time-Reversal  Array  which  can  be  used  for  the  passive  iterative  time- 
reversal  process. 

Only  a  rough  estimate  of  the  trajectory  is  necessary  to  estimate  the  relative 
position  of  the  synthetic  aperture  TRA  elements  with  respect  to  the  target  location.  This 
allows  one  to  estimate  which  bottom  patch  corresponds  to  a  specific  given  time-gated 
window  of  the  backscatter  signals  in  order  to  perform  a  search  for  the  target  along  the 
seafloor  (see  Fig.  4).  Reflectivity  Maps  (similar  to  Fig.  2)  over  the  whole  search  grid 
could  then  be  generated  to  help  detect  and  localize  the  target. 

Validation  of  the  detection  capabilities  of  the  synthetic  aperture  time-reversal 
processing  in  this  simple  bi-static  geometry  would  provide  an  initial  demonstration  for 
self-adaptive  and  robust  MCM  application  eventually  involving  AUVs.  The  data/analysis 
done  in  the  near  field  by  Schmidt  et  al  will  also  be  studied  in  the  context  of  relating  the 
decomposition  of  the  TR  process  to  the  complex  elastic-properties  of  the  combination  of 
target- sediment. 
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Figure  4:  Top  view  implementation  for  the  Target  using  a  synthetic  aperture  Time- 
Reversal  Array.  Search  grid  along  the  the  seafloor  can  be  constructed  using  an  estimated 
position  of  the  Synthetic  aperture  TRA  with  respect  to  the  ocean  bottom,  based  on  simple 
navigation  data.  Reflectivity  map  similar  to  Fig.  2  can  then  be  constructed  over  the  whole 
search  grid. 


SWAMSI-2nd  Experiment 

Although  beyond  the  time  frame  of  this  proposal,  it  is  likely  that  a  second  SWAMSI 
experiment  will  be  conducted  in  FY08.  The  algorithms  developed  here  will  be  available 
for  use  in  the  FY08  experiment  which  will  involve  multiple  platforms. 


Personnel:  In  addition  to  the  Pi’s,  there  will  be  some  programming/data  work  provided 
by  the  appropriate  MPF  staff. 


Results:  The  results  of  this  project  have  led  to  new  algorithms  for  SWAMSI  experiments 
and  are  summerarized  in  the  follow  extract  from  W.  Higley’s  thesis: 


Chapter  from  W.  Higley's  Thesis 


Chapter  5 

Detection  of  Resonances  on  a  Spherical  Target 

abstract 

The  use  of  low  frequency  (2-15  kHz)  acoustic  signals  for  mine  eountermeasures  have 
been  motivated  by  the  ehance  to  use  scattering  features  of  buried  targets  (e.g.  the  radius, 
wave  thiekness,  and  elastie  properties)  in  the  elassifieation  of  the  target.  At  frequeneies, 
lower  than  traditionally  used,  sound  penetrates  into  the  sediment  better  and  excites  less 
signal-generated  noise  off  the  seabed.  The  struetural  waves  of  typieal  mine-sized  targets 
are  enhaneed  at  these  frequencies.  This  paper  introduces  and  compares  two  forms  of 
proeessing  designed  to  elassify  a  target  as  resonant  or  non-resonant  based  on  the 
baekseattered  field.  The  two  forms  of  proeessing  are  distinguished  as  an  ineoherent, 
energy-based  method,  and  a  eoherent  method  based  on  the  eoneept  of  prineipal 
eomponent  analysis  (PC A). 

5.1  Introduction 

When  performing  mine  eountermeasures  (MCM),  it  is  preferable  to  have  a  large 
standoff  distanee.  That  is,  it  is  desirable  to  deteet  and  elassify  a  mine  while  as  far  away 
as  possible  from  it.  This  desire  leads  to  the  use  of  lower  frequeneies,  whieh  propagate 
further  horizontally.  This  decreases  the  resolution  of  traditional  aeoustic  imaging 
teehniques.  Also,  the  inerease  in  standoff  distance  leads  to  an  inerease  in  the  infiuenee  of 
the  waveguide  in  introducing  multipath. 


103 


Chapter  from  W.  Higley's  Thesis 


104 


Recent  experiments  [1],[2]  have  attempted  to  use  low  frequency  (2-15  kHz)  signals  as 
part  of  the  target  detection  and  classification  problem  of  MCM.  At  this  low  frequency, 
although  the  resolution  of  traditional  imaging  is  limited,  other  properties  such  as  the 
radius,  wall  thickness,  and  elastic  properties  of  the  target  can  be  extracted. 

Lower  frequencies  also  benefit  from  greater  seabed  penetration.  Attenuation  above 
the  grazing  angle  decreases  approximately  linearly  with  frequency,  and  below  the  grazing 
angle  the  evanescent  lateral  wave  field  in  the  bottom  penetrates  deeper  [3].  The  signal¬ 
generated  noise  from  reverberation  off  surface  features  of  the  seafloor  (such  as  sand 
ripples)  will  also  decrease  with  frequency  [4]. 

The  problem  of  using  target  resonances  to  classify  a  target  is  complicated  in  a 
waveguide  by  the  multipath  in  the  environment.  This  paper  compares  two  forms  of 
processing  on  two  targets  in  both  freespace  and  in  a  waveguide.  The  two  forms  of 
processing  are  classified  as  an  incoherent,  energy-based  processing  method,  and  a 
coherent  method  based  on  the  concept  of  principal  component  analysis  (PCA).  Our 
simulations  show  the  benefit  of  the  coherent  method  over  the  incoherent  method  in 
classifying  a  target  as  resonant  or  non-resonant.  Finally,  this  paper  proposes  an 
experiment  using  a  horizontally-aligned  array  towed  from  an  autonomous  underwater 
vehicle  (AUV). 

5,2  Freespace  Simulation 

A  freespace  simulation  is  done  to  test  the  feasibility  of  our  processing.  In  an  aqueous 
medium  with  a  density  of  1  g/cm  and  a  sound  speed  of  1507  m/s,  there  is  an  array  of 
nine  hydrophones  arranged  vertically  with  a  separation  distance  of  10  m.  The  target  is 
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Figure  5.1:  The  geometry  of  the  simulation.  A  vertieal  array  eonsisting  of  nine  elements 
is  oriented  vertieally  in  freespaee.  The  distanee  between  eaeh  element  is  10  m,  and  the 
first  element  is  eonsidered  10  m  below  some  arbitrary  origin.  The  target  is  loeated  sueh 
that  the  bottom  of  the  target  is  10  m  below  the  lowest  element,  and  200  m  away  from  the 
vertieal  array. 

situated  sueh  that  the  bottom  is  loeated  10  m  below  the  lowest  souree/reeeiver,  and  the 
eenter  is  loeated  200  m  away  horizontally.  The  geometry  of  the  set-up  is  shown  in  Fig. 
5.1.  Two  targets  are  analyzed:  a  bubble  and  a  spherieal  shell.  The  bubble  is  a  spherieal 
vaeuum  with  a  radius  of  .53  m.  The  spherieal  shell  has  an  inner  radius  of  .5  m  and  an 
outer  radius  of  .53  m.  The  interior  of  the  shell  is  a  vaeuum,  and  the  shell  itself  is 
eomposed  of  steel  with  a  density  of  7.7  g/em  ,  a  shear  sound  speed  of  2940  m/s  and  a 
eompressional  sound  speed  of  5880  m/s.  The  two  targets  are  shown  in  Fig.  5.2. 

Eaeh  hyrdrophone  transmits  and  pulse  whieh  is  baek-seattered  by  the  target  and 
reeeived  on  all  nine  hyrdrophones.  This  ereates  a  three-dimensional  matrix  for  eaeh 
target  as  a  funetion  of  transmitter  depth,  souree  depth,  and  frequeney.  The  matriees  for 
the  bubble  and  the  sphere  are  denoted  as  Kh(zt,Zs,co)  and  Ks(zt,Zs,co),  respeetively.  By 
taking  the  inverse  Fourier  transform  along  the  frequeney  domain,  one  ean  examine  the 
reeeived  field  for  any  given  transmitter.  In  Fig.  5.3,  the  freespaee  responses  from  the  first 
transmitter  for  both  the  bubble  and  the  shell  are  shown.  The  response  from  the  bubble  is 
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p  =  7.7  g/cnri^ 


Figure  5.2:  The  two  targets,  both  loeated  in  an  aqueous  medium.  The  first  target  is  a 
bubble,  whose  interior  is  a  vaeuum.  The  radius  of  .53  m  is  matehed  to  the  outer  radius  of 
the  seeond  target,  a  resonant  shell.  The  wall  thiekness  of  the  shell  is  3  em,  and  it  is  also 
filled  with  a  vaeuum.  The  speed  of  sound  and  density  of  the  water  are  1507  m/s  and  1 
g/em^  respeetively.  The  shell  is  eomposed  of  steel  with  a  density  of  7.7  g/em^.  The 
sound  speeds  of  eompression  and  shear  waves  in  the  steel  are  5880  m/s  and  2940  m/s 
respeetively. 


composed  of  only  a  specular  component,  that  field  that  reflects  off  the  surface  of  the 
bubble.  Flowever,  the  response  from  the  shell  is  composed  of  the  same  specular 
component  (possibly  phase-shifted)  and  an  additional  resonant  component  that  arrives 
later. 

Given  the  matrices  lLb{zt,Zs,co)  and  Ks(z/,Zi,m),  it  is  desirable  to  find  some  processing 
that  will  distinguish  between  the  two.  As  can  be  seen  in  Fig.  5.3,  in  the  freespace  case, 
one  could  simply  time-gate  out  the  specular  and  be  left  with  the  resonance,  if  it  existed. 
However,  this  method  will  not  work  in  a  waveguide,  as  will  be  seen  later,  so  other 
processing  methods  need  be  examined. 

The  first  method  considered  is  an  incoherent,  energy-based  method.  The  matrix  K 
can  be  reduced  to  a  one-dimensional  vector  by  summing  incoherently  across  transmitters 
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Figure  5.3:  The  baekseattered  response  on  the  vertieal  array  from  both  the  bubble  (gray) 
and  the  shell  (blaek)  when  exeited  by  a  single  souree  on  the  array.  The  response  from  the 
bubble  eonsists  of  a  strong,  single,  speeular  eeho,  while  the  response  from  the  shell  is 
eomposed  of  both  that  same  speeular  response,  and  a  resonant  eomponent  that  arrives 
later. 


and  reeeivers.  The  proeessing  ean  be  expressed  mathematieally,  where  N  is  the  number 
of  hydrophones,  as 

N  N 

z^=l  z,=l 

This  term  ean  also  be  expressed  in  terms  of  the  singular  value  deeomposition  (SVD)  [5]. 
Eaeh  frequeney  of  K  represents  an  NxN  matrix.  Taking  the  singular  value  deeomposition 
of  this  matrix  yields  N  singular  values.  Thus,  there  are  N  singular  values  at  eaeh 
frequeney.  The  sum  of  these  singular  values  squared  also  yields  the  veetor  e(co).  This  is 
shown  mathematieally  below  as  [6]: 

svd  (K(z,,z^,  0}))  ^  cr.  (cd)  for  /=  1 ,2, . . .  ,A^ 

N  N  N 

X  =  Z  Z  • 

i=l  Z^=lZj=l 


(5.2) 
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Figure  5.4:  The  reeeived  energy  aeross  frequeney  of  the  baekseattered  response  of  the 
bubble  (gray)  and  the  shell  (blaek),  averaged  aeross  souree  and  reeeiver.  A  few  of  the 
resonant  peaks  of  the  shell  are  seen  at  lower  frequeneies. 

Shown  in  Fig.  5.4  are  the  results  of  this  energy-based  proeessing.  The  resonanees  of 
the  shell  are  seen  as  peaks  in  the  average  energy  at  speeifie  frequeneies.  The  bubble 
shows  no  sueh  peaks  in  the  frequeney  domain,  and  it  rather  flat  aeross  frequeney,  as  one 
would  expeet  from  a  speeular  eeho. 

The  seeond  method  eonsidered  is  a  eoherent  one  related  to  prineipal  eomponent 
analysis.  We  reduee  the  three-dimensional  matrix  K(z^,z^,(u)  into  a  two-dimensional 

matrix,  labeled  K(z^,z^  -co) ,  via  eoneatenation.  The  reeeived  field  in  the  frequeney 

domain  from  eaeh  souree  ean  be  thought  of  as  a  separate  observation  of  a  random  veetor. 
From  these  observations,  a  sample  eovarianee  matrix  ean  be  ereated  as  follows: 

1  ^  , 

"  TTT  ^  ^)  • 

^  z,=l 


(S.3) 
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Taking  the  first  eigenvector  of  the  matrix  S  yields  the  first  principal  component  of  the 
data  [6].  The  first  principal  component  is  the  energy-constrained  vector  which 
maximizes  the  average  correlation  with  the  data.  It  can  be  thought  of  as  the  part  of  the 
data  that  is  the  most  similar  across  all  observations,  or  in  the  case  examined  here,  across 
all  sources.  We  write  the  first  principal  component  mathematically  as; 


p(z^  •  0))  =  arg  max  ^  p"  (z^  •  ni)K(z, ,  z^  •  ni) 


z,=l 


(S.4) 


As  stated  above,  p  is  also  the  first  eigenvector  of  the  sample  covariance  matrix,  S.  The 
vector  p  can  be  found  another  way,  using  singular  value  decomposition.  If  we  consider  S 
to  have  the  following  eigenvalue  decomposition, 

S  =  UAU",  (S.5) 

we  can  write  the  singular  value  decomposition  of  K  as 

K  =  UV(A^-1)AV".  (S.6) 

In  both  cases,  the  matrix  U  is  identical.  The  principal  component,  p,  can  be  thought  of  a 
either  the  first  eigenvector  of  S,  or  the  first  left  singular  vector  of  K  . 

Deconcatenation  of  the  vector  p(z^  •  co)  results  in  a  matrix  P(zr,co)  that  is  a  function  of 

receiver  and  frequency.  Taking  the  inverse  Fourier  transform  along  the  frequency 
domain  yields  a  single  time-domain  field  for  each  matrix,  K/,  and  K,.  These  fields  are 
plotted  in  Fig.  5.5.  The  fields  shown  below  can  be  thought  of  as  the  parts  of  the  received 
fields  that  are  most  independent  of  source  depth.  In  the  case  of  the  bubble,  there  is  a 
strong  first  arrival  as  there  is  similarity  between  the  specular  echoes  for  each  source.  In 
the  case  of  the  shell,  the  resonance  is  nearly  separated  from  the  specular,  as  the  resonance 
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Figure  5.5:  The  first  prineipal  eomponent  from  the  bubble  (gray)  and  the  shell  (blaek). 
The  eoherent  proeessing  first  yields  a  field,  whieh  is  shown  in  these  figures.  The  field 
extraeted  from  the  bubble  is  eomposed  of  a  strong  first  arrival  and  weaker  later  arrivals 
whieh  may  be  attributed  to  weak  resonanees  of  the  bubble.  On  the  other  hand,  the 
extraeted  field  from  the  shell’s  response  is  eomposed  of  periodie  pulses  eorresponding  to 
the  resonanees  propagating  around  the  shell. 


is  entirely  souree  independent.  The  resonant  eomponent  of  the  reeeived  field  is 
independent  of  exeitation,  and  shows  up  prominently  is  this  proeessing. 
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Figure  5.6:  The  result  of  the  eoherent  proeessing  on  the  bubble  (gray)  and  the  shell 
(blaek).  The  periodie  resonanee  peaks  of  the  shell  are  quite  pronouneed  through  the 
frequeney  range,  whereas  the  response  of  the  bubble  is  relatively  flat,  as  one  would 
expeet  for  a  speeular  eomponent. 


We  ean  reduee  the  matrix  P(z^,co)  into  a  frequeney  veetor  by  summing  ineoherent 
aeross  reeeiver.  We  ean  eompare  this  new  veetor,  q(co),  with  the  previously  ealeulated 
ineoherent  veetor,  The  new  veetor  is  written  mathematieally  as: 

q(n>)  =  f;|P(z,,^y)|^  (S.7) 

z^=l 

The  veetor  q(co)  is  ealeulated  for  both  the  bubble  and  the  sphere,  with  the  results  shown 
in  Fig.  5.6.  The  proeessing  requires  the  results  for  both  the  bubble  and  the  sphere  to  have 
the  same  magnitude,  when  summed  aeross  frequeney,  as  the  matrix  U  is  unitary, 
regardless  of  the  relative  energy  in  the  matrix  K.  The  result  of  the  prineipal  eomponent 
analysis  proeessing  done  on  the  bubble  is  relatively  flat,  eompared  to  the  results  obtained 
when  done  on  the  shell.  The  resonanees  of  the  shell  appear  again  as  peaks  in  the 
frequeney  domain,  with  a  spaeing  of  approximately  675  Hz.  This  eorresponds  to  a 
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Figure  5.7:  The  result  of  baek-propagating  the  field  extraeted  in  Fig  5.5b.  The  foeus  is 
seen  only  at  the  loeation  of  the  target,  at  a  depth  of  100  m  and  a  range  of  200  m. 

repeated  time-delay  of  about  1.5  ms,  a  spacing  that  is  seen  in  Fig  5b.  This  is  the  time  it 
takes  the  resonance  to  travel  around  the  sphere  before  emitting  again.  From  this,  we  can 
calculate  the  group  speed  of  the  Aq  wave  as  approximately  2120  m/s. 

Comparing  Fig.  5.4  and  Fig.  5.6,  we  notice  that  the  principal  component  analysis  is 
able  to  distinguish  resonant  peaks  at  10  kHz  and  above,  whereas  the  energy-based 
analysis  cannot.  Also,  those  peaks  at  frequencies  less  than  10  kHz  are  more  pronounced 
than  in  the  energy-based  analysis. 

Returning  our  attention  to  the  time-domain  resonance  field  of  Fig.  5.5b,  expressed  by 
the  matrix  Vs{zr,co),  an  incoherent  back-propagation  of  this  field  can  be  calculated.  The 
result  of  this  calculation  is  shown  in  Fig.  5.7.  We  see  that  the  back-propagated  resonance 
field  focuses  at  a  location  just  above  the  sphere.  The  relative  tightness  of  the  focus  is 
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Figure  5.8:  a)  In  gray  is  the  field  previously  extraeted  from  the  shell  and  shown  in  Fig. 
5.5b.  The  windowed  portion,  shown  in  blaek,  eorresponds  to  a  resonant  emission  from 
the  target,  b)  The  result  of  baek-propagating  the  resonant  emission  shows  a  pair  of  foei 
both  above  and  below  the  loeation  of  the  target,  eorresponding  to  the  two  loeations  on  the 
target  from  whieh  resonanees  are  emitted.  The  other  two  foei,  in  front  of  and  behind  the 
target,  are  assumed  sidelobes,  artifaets  of  the  proeessing. 


related  to  the  vertieal  aperture  of  the  hydrophone  array  used.  Baek-propagation  of  the 
matrix  P/,(z^,(u)  results  is  a  similar  foeus  near  the  loeation  of  the  sphere.  A  differenee 
between  the  two  baek-propagated  fields  ean  be  found  by  windowing  out  a  single 
resonanee  of  the  matrix  P^(z^,co)  in  the  time  domain  (as  shown  in  Fig.  5.8a)  and  baek- 
propagating  only  that  pieee.  As  seen  in  Fig.  5.8b,  this  results  in  a  group  of  foei  in  the 
area  around  the  sphere.  The  two  middle  foei  are  interpreted  at  the  two  patehes  on  the 
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sphere  from  whieh  resonanees  are  emitted  to  the  array.  The  other  two,  in  front  and 
behind  the  sphere,  are  interpreted  as  sidelobes,  artifaets  of  the  proeessing. 

1,3  Waveguide  Simulation 

Having  established  two  forms  of  proeessing,  an  ineoherent  energy-based  method 
yielding  e(co),  and  a  eoherent  method  based  on  prineipal  eomponent  analysis  whieh 
yields  q(co),  the  ease  of  a  buried  sphere  in  a  waveguide  is  now  examined.  The  array  is 
situated  vertieally  as  before,  in  a  100-m  deep  underwater  waveguide.  The  target  is 
situated  200  m  away  horizontally,  as  before,  but  instead  of  the  bottom  of  the  target  being 
at  a  depth  of  100  m,  the  top  is  loeated  at  that  same  position.  The  target  is  therefore 
buried,  sueh  that  the  top  of  the  target  is  tangent  to  the  waveguide  bottom.  The  noise-free 
matriees  Kft(Z(,Zi,co)  and  K,(Z(,Zi,co)  are  aequired  as  before,  and  an  additional  matrix 
Kr(zt,Zs,co)  is  ealeulated  from  a  seafloor  reverb  model,  to  allow  us  to  examine  our 
proeessing  in  the  presenee  of  signal-generated  noise.  The  baekseattered  response  of  the 
flush-buried  sphere  in  a  waveguide  is  shown  in  Fig.  5.9  for  the  two  targets;  the  bubble 
and  the  sphere.  The  shape  of  the  baek-seattered  response  is  the  same,  but  eaeh  pulse  of 
the  bubble’s  response  is  eonvolved  with  the  resonanee  in  the  ease  of  a  shell. 

The  two  above-mentioned  method  of  proeessing  are  applied  now  to  the  ease  of  a 
flush-buried  target.  For  eomparison  purposes,  eaeh  of  the  veetors  e(co)  and  q(co)  are 
normalized  sueh  that  the  energy  in  eaeh  veetor  is  one,  whieh  we  denote  e(co)  and  q{co) . 
Figure  5.10  shows  the  result  of  both  forms  of  proeessing,  without  noise,  on  both  flush- 
buried  targets.  In  the  ease  of  the  bubble  (Fig.  5.10a),  there  are  no  notieeable  peaks  found 
in  either  form  of  proeessing,  as  we  would  expeet  from  speeular  eehoes.  Both  the 
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Figure  5.9:  The  baek-seattered  response  from  a  flush-buried  bubble  (gray)  or  shell  (blaek) 
in  a  waveguide.  The  shape  of  the  response  in  the  same  in  both  eases,  but  eaeh  pulse  of 
the  response  in  the  ease  of  the  bubble  is  turned  into  a  resonant  arrival  strueture  in  the  ease 
of  the  shell. 


ineoherent  and  eoherent  forms  of  proeessing  yield  resonant  peaks  when  applied  to  the 
flush-buried  shell  (Fig.  5.10b),  but  the  method  based  on  prineipal  eomponent  analysis 
eauses  these  peaks  to  be  more  pronouneed. 
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Figure  5.10;  a)  The  result  of  the  ineoherent  (gray)  and  eoherent  (blaek)  proeessing  on  the 
non-resonant  bubble  in  a  waveguide.  The  frequeney  response  is  relatively  flat,  as  we 
would  expeet  from  a  purely  speeular  eeho.  b)  The  result  of  the  ineoherent  (gray)  and 
eoherent  (blaek)  proeessing  on  the  resonant  shell  in  a  waveguide.  The  periodie  resonant 
peaks  are  elearly  identified,  and  eorrespond  well  to  those  found  in  freespaee  (Fig.  5.6). 


By  examining  the  results  of  the  proeessing  under  various  levels  of  signal-generated 
noise,  a  differenee  is  again  seen  in  the  results  of  the  ineoherent  and  eoherent  proeessing. 
Figure  5.11  shows  the  results  of  both  forms  of  proeessing  on  the  flush-buried  shell  for 
four  different  noise  levels.  In  Fig.  5.11a,  the  results  are  very  similar  to  Fig.  10,  with  both 
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Figure  5.1 1:  A  comparison  of  the  incoherent  (gray)  and  coherent  (black)  processing  on 
the  flush-buried  shell  for  four  different  signal-to-noise  ratios.  Although  both  methods 
degrade  as  noise  increases,  the  coherent  method  still  shows  peaks  at  -10  dB,  whereas  the 
incoherent  method  fails  at  this  SNR. 


forms  of  processing  showing  strong  resonant  peaks,  and  those  of  the  coherent  PCA-based 
method  a  little  stronger.  In  Fig.  5.11b  and  5.11c,  a  degradation  of  performance  for  both 
methods  is  witnessed,  with  the  coherent  method  still  performing  better  than  the 
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incoherent  energy-based  method.  In  Fig.  5.1  Id,  the  ineoherent  method  fails,  whereas 
small  peaks  at  the  resonant  frequeneies  are  still  seen  using  the  eoherent  method. 

The  main  eonelusion  drawn  from  Fig.  5.10  and  5.11  is  the  better,  more  peaked 
performanee  of  the  eoherent  method,  q((y) ,  over  the  ineoherent  method,  e((y) .  These 
results  ean  be  distilled  into  a  elassifieation  metrie  by  ealeulating  the  sample  varianee  or 
sample  kurtosis  of  the  above  veetors  and  eomparing  these  values  at  various  signal-to- 
noise  ratios  to  the  varianee  or  kurtosis  of  the  same  proeessing  performed  on  a  bubble. 
Figure  5.12  shows  the  result  of  just  sueh  a  proeessing.  Although  the  sample  varianee  the 
eoherent  proeessing  (blaek  in  Fig.  5.12a)  is,  in  general,  higher  than  that  of  the  ineoherent 
proeessing  (gray  in  Fig.  5.12a)  the  separation  of  the  eoherent  method  when  applied  to  the 
shell  (solid  line)  and  the  bubble  (dashed  line)  is  greater  than  the  separation  of  the 
ineoherent  method  applied  to  the  shell  and  the  bubble.  The  eoherent  method  reeognizes  a 
differenee  between  the  two  down  to  an  SNR  of  -10  dB,  whereas  the  ineoherent  method 
fails  at  this  reverb  level.  This  result  was  previously  seen  in  Fig.  5.1  Id.  Sinee  the 
eoherent  PCA-based  method  yields  a  more  peaked  result,  the  sample  kurtosis  of  e{a>) 
and  q((y),  shown  in  Fig.  5.12b,  yields  another  elassifieation  metrie.  In  this  ease,  the 
eoherent  method  reeognizes  a  differenee  between  the  shell  and  the  bubble  at  SNRs  as  low 
as  -20  dB,  whereas  this  metrie  is  rather  poor  for  the  ineoherent  method,  failing  at  an  SNR 


of -5  dB. 
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5.4  Experimental  Proposal 

Although  both  forms  of  processing  perform  well  when  a  vertieal  array  is  used  at  high 
signal-to-noise  levels,  sueh  a  eonfiguration  would  be  impraetieal  to  realize  at-sea  for  real- 
world  applieation.  Of  great  interest  would  be  the  results  of  this  proeessing  performed  on 
a  dataset  realized  from  a  horizontal  array  towed  from  an  AUV. 

The  strength  of  the  proeessing,  in  both  the  ineoherent  and  eoherent  ease,  is  based  on 
“averaging”  away  the  frequeney-domain  peaks  eaused  by  waveguide  interferenee, 
leaving  only  those  peaks  emitted  from  a  resonant  target.  Therefore,  sueh  a  towed  souree 
would  have  to  sample  the  wavenumber  spaee  effieiently,  whieh  is  done  mueh  more  easily 
with  a  vertieal  array  than  a  horizontal  one  [7].  The  best  solution  to  this  would  seem  to  be 
a  broadside  array  as  elose  as  feasible  to  the  target.  The  point  of  aeoustie  deteetion  and 
olassifieation  of  a  target  in  MCM,  of  eourse,  is  to  not  be  elose  to  the  target,  so  a 
eompromise  in  distanee  must  be  made. 

5.5  Conclusion 

This  paper  has  introdueed  and  eompared  two  forms  of  proeessing  to  elassify  a  target 
as  either  resonant  or  non-resonant  by  using  the  baek-seattered  field  from  that  target.  The 
two  proeessing  methods  are  an  ineoherent,  energy-based  method,  and  a  eoherent  PCA- 
based  method.  Comparison  of  the  two  methods  in  simulation  shows  the  benefit  of  the 
eoherent  method  over  the  ineoherent  method  in  deteeting  resonanees  in  the  presenee  of 
noise.  Lastly,  an  experiment  to  test  the  praetieality  of  the  proeessing  for  real-world 
MCM  is  proposed. 
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Figure  5.12:  a)  Sample  variance  of  the  frequency  vectors  corresponding  to  the  incoherent 
processing  (gray)  and  coherent  processing  (black).  The  solid  lines  correspond  to  the 
processing  performed  on  the  flush-buried  shell,  and  the  dashed  lines  correspond  to  the 
processing  performed  on  the  bubble  for  comparison.  The  coherent  processing 
outperforms  the  incoherent  processing  by  approximately  5  dB,  using  this  metric,  b)  The 
sample  kurtosis,  which  is  high  for  peaky  vectors,  of  the  frequency  vectors  corresponding 
to  the  incoherent  processing  (gray)  and  coherent  processing  (black).  The  solid  lines 
correspond  to  the  processing  performed  on  the  flush-buried  shell,  and  the  dashed  lines 
correspond  to  the  processing  performed  on  the  bubble  for  comparison.  Using  this  metric, 
classification  using  the  coherent  PCA-based  method  is  possible  as  low  as  -20  dB. 
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